The density functional theory of classical fluids revisited. 
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We reconsider the density functional theory of nonuniform classical fluids from the point of view 
of convex analysis. From the observation that the logarithm of the grand-partition function logH[(/3] 
is a convex functional of the external potential tp it is shown that the Kohn-Sham free energy A[p\ is 

■ a convex functional of the density p. logE[(/p] and A[p\ constitute a pair of Legendre transforms and 
O ' each of these functionals can therefore be obtained as the solution of a variational principle. The 

convexity ensures the unicity of the solution in both cases. The variational principle which gives 
logH[(/j] as the maximum of a functional of p is precisely that considered in the density functional 
theory while the dual principle, which gives A[p\ as the maximum of a functional of seems to be 
, a new result. 

' KEY WORDS: Density fonctional theory; Theory of liquids; Convex analysis. 

B 

^ ■ I. INTRODUCTION 

o 

I ^ I The convexity properties of the thermodynamic potentials play an important role in thermodynamics. For instance, 

in the case of a simple fluid, the entropy S{U, V, N) is a concave function of the internal energy U, the volume V 
and the number of particles N. As discussed by Wightmanld this property, supplemented by the hypothesis that the 
— ^ ^ entropy is a homogeneous function of {U, V, N) of the first degree (i. e. S{XU,^V, XN) — XS{U, V, N) for any A > 0), is 
fully equivalent to the second law of thermodynamics in Gibbs' formulation.0 However, it should be stressed that the 
concavity of S" as a function of {U, V, N) - and, similarly, the convexity of C/ as a function of {S, V, AQ- are properties 
, which hold only in the thermodynamic limit where the properties of extensivity of S and U are valid.cl However, there 

■ is one family of convexity properties which holds even in a finite volume. That is the convexity of the logarithm of 
, the canonical partition function log 2' as a function of the inverse temperature f3 — \/kT, and the convexity of the 

■ logarithm of the grand partition function log S in /? and v = (/x chemical potential) a 
In this paper we extend the latter property to the case of nonuniform systems. The system under interest is thus a 

classical fluid in a volume V and in the presence of an external one-body potential (p{x). It will be shown in Sec. HI 
that log S is a convex functional of the generalized potential 



> 



u{x)^iy-P^ix). (1.1) 

O ' . . . . 

: The implications of this mathematical property are studied in details and are shown to lead to a new formulation of 
the density functional theory (DPT) of classical fluids. Indeed, the convexity of log S [u] implies the existence and the 
unicity of a convex functional f3A [p] of the inhomogeneous density p (x) such that log S [u] and (3 A [p] are a couple 
of Legendre transforms. Of course, the functional (3A [p] identifies with the intrinsic Helmholtz (Kohn-Sham) free 
5h [ energy. Since they are Legendre transforms, both log S [u] and jSA [p] can be expressed as the solution of a variational 
principle. The variational principlerwhich defines log S [u] as the maximum of a fuactjonal of p{x) is precisely the one 
which was noticed by Kohn et aLSa and constitutes the cornerstone of the DFTHnZl The dual variational principle, 
which defines (3A [p] for a given p (x) as the maximum of a functional of u{x) seems to have been overlooked. 
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Our paper is organized as follows. In Sec. II we review some mathematical properties of the grand-partition function 
of a nonuniform fluid. In Sec. Ill we show that logS [u] is a convex functional of the generalized potential u{x) which 
allows us to define its Legendre transform A [p] . The link with the DFT is discussed in details. Conclusions are drawn 
in Sec. IV. We have found usefull to give easy (formal) proofs of the theorems of convex analysis which are needed in 
the text in an appendix. 



II. THE GRAND PARTITION FUNCTION AND ITS FUNCTIONAL DERIVATIVES 



In this section we review known properties of nonuniform classical fluids and fix our notations. We consider a 
simple fluid made of identical, structureless classical particles in a finite volume V. The dimension of space is d. A 
configuration of the system in the grand-canonical ensemble will be denoted by w = {N, xi, . . . -Xn) where Xi denotes 
the position of the ith particle. The grand-canonical partition function S [u] can be written asU 



E[u]^ I dA*Hexp((/5H|u)) , (2.1) 
where we have introduced the scalar product 

{p{uj)\u) = I d'^x p{x;u})u{x) (2.2) 



between the generalized potential u{x) (cf Eq. (1.1)) and and the microscopic density of particles 



N 



p{x; w) = ^ 5'^[x - Xi) 



i=l 



In Eq. (2A) we have adopted the compact notation 

dniuj) = duj exp {-pW{u;)) , (2.3) 

where W{uj) denotes the internal potential energy of N particles in the configuration lu and dio is the grand-canonical 
measure defined as 

oo 

'^^^Ejmj^ d'^i---d'xr,. (2.4) 

where A is De Broglie thermal wave length. The grand-canonical average of any microscopic variable A{uj) can thus 
be written as 

{A)= j dp{uj)eyi^{{p{uj)\u))A{uj)/ ^[u] . (2.5) 

We shall denote asgusual by U,[u] = — /3^^1ogS[u] the grand Potential; note that it can be rewritten as a grand- 
canonical averageOi3 

(}n[u] = {pW{u;) + log /o(co)) - {p\u) , (2.6) 

p{x) ^ {p{x-u)) , (2.7) 

where /o(w) = exp {—l3W{uj) + {p{uj)\u)) is the grand-canonical equilibrium distribution in the phase space. 

We shall denote by U the set of potentials u{x) such that the r.h.s. of Eq. (2.1) converges. For the sake of simplicity 
we shall restrict ourselves to the case of i?-stable systems in the sense of Ruelle, i.e. systems such that W{u}) > —NB 
where the lower bound B < oo \s some constant independent of N . Most models of interest satisfy this property. We 
have obviously 

S[u] < exp ( -^exp(/3S) / d''xexp(M(f)) ) , (2.8) 



met 

system; i.e. W = {u : V ^ M ; exp(u) g Ly [d'^a?]} . It can be shown that, for H- stable systems, W is a convex set 



from which it follows that U is the set of the functions u{x) such that exp {u(x)) is integrable over the volume V of thj8| 

" ~ tl 
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Indeed, let uq and ui two potentials belonging to U. We define the closed (respectively open) interval [ito,wi] as the 
set of functions u\ = {1 — X)uo + Xui where < A < 1 (respectively < A < 1). As a consequence of the convexity 
of the exponential function exp ((1 — X)uo + Xui) < (1 — A) exp(Mo) + Aexp(Mi) and therefore [liiftrffJ C U, QED. 

As it is well known, S[u] is the generating functional of a hierarchy of correlation functions.uilJ'LlI Expanding the 
expression (2.1) of S[w + Su] in powers of Su one finds for {u, u + Su € U) 

°° 1 f 

E[u + Su]/ E[u] = V -r / d'^xi ■ ■ ■ d'^Xn G'("Hfi, . . . , f„) Su{xi) ■ ■ ■ Ju(f„) , (2.9) 

n=0 

where the correlation function G^"^ has been defined as 

1 E[u] 



G(")(fi,...,af„) = - 



^[u] Su{xi) • • • 5u{Xn) 

71 

= {l[p{xf,cj)) . (2.10) 

1=1 

It is sometimes also useful to introduce the connected correlation functions as 

gr^(x....,a-.o- ,f:V°fr^ - (2.11) 

OU(Xi) ■ ■ ■ du[Xn) 
The relation between ci"'' and G^"^ can be written symbolically 

Gl")(fi,...,f„)=G(")(fi,...,f„)-^ n Gi'^Hx,„...,x,J , (2.12) 

m<n 

where the sum of products is carried out over all possible partitions of the set {!,... ,n) into subsets of cardinal 
m < n.E3EJ For instance we have 

G(2)(fi,f2) =G(2)(a-i,f2)-p(fi)p(f2) , (2.13) 

where p = G^^^ = ci^^ is the mean density of particles. Note that in the theory of liquids one rather defines the 
correlations functions p'"-* as the functional derivatives of S with respect to the activity z — exp (71). More precisely : 

^ ) - nLi^(^-o ^(")iogs[z] 

^ (^i"--'^"^- E[z] 5z{x,)---5z{x^)- ^'-'^^ 
The functions and G^"-* differ by delta functionsil0, for instance G'^''(a;i,X2) = P^'^\xi,X2) + p(xi)5'^{xi — X2). 

III. CONVEXITY AND LEGENDRE TRANSFORMATIONS 

A. Convexity of logH[u] 

Let us show now that logS[u] is a strictly convex functional of u, i.e. that, for any ui and U2 ^ ui E U we have, 
for any < A < 1 the strict inequality 

logS[AMi + (1 - A)w2] < AlogS[ui] + (1 - A)logSH ■ (3.1) 

Let u{x) be an arbitrary potential of U and h{x) any non-zero function defined on the volume of the system. The 
quadratic form 



(logS)^^^ [h] ^ / d''x,d^x2 . ,^ " ,l\ Kx.Mx,) (3.2) 
Jv du{xi)du[x2) 

will be shown to be str ictly positi ve which implies the strict convexity from theorem II of the appendix. Indeed it 



follows from Eqs. (2.10) and (2.11) that 
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{\ogE)\^> [h] = / d''xid''x2 Gl^\xi,X2)h{xi)h{x2) , 



(H^) - {Hf > 



(3.3) 



where H{lo) = Y^iLi H^i)- 

One can also give an alternative, more direct proof of the convexity of logS[w]. For ui, M2 ^ wi and < A < 1 
we have, with slightly simplified notations, 



E[Xui + (1 - A)u2] 



d^i exp(A(/5|ui))exp((l - A)(p|u2)) 



< 



dfi [exp(A(/5|ui))]^/^ 



1 ^ 




X 


[X 







(l-A) 



(3.4) 



which yields inequality (3.1) by taking the logarithm. In order to write the second line of Eq. (3.4) we have made-ase 
of Holder's inequality which states that, for any real numbers p, q such that p > 1, <? > 1, l/p + 1/q = 1 one haai3 



dfiiuj) f{uj)g{uj)\ < / dpi{Lo)\f{uj)\P 



i/p 



dfi{uj)\g{uj)\' 



1/9 



(3.5) 



The case of equality in Eq. ( |3.5[ ) occurs only if |/(w)|p/|(7(w)|'^ is almost everywhere equal to a constant. This remark 
implies the strict inequality in Eq. ( |3.1[ ). Note that the grand-potential /3ri[u] = — logS[u] is therefore a strictly 
concave functional of u{x) even before the passage to the thermodynamic limit (TL). 



B. The Kohn-Sham free energy as the Legendre transform of log 2 



Let us define the functional 



PA[p,u] = {p\u) - log^[u] 



(3.6) 



For a given function p{x), l3A[p, u] is obviously a concave functional of u. Consequently, if f3A[p, u] has a maximum 
for some u £ U this maximum is strict and therefore unique. This is theorem III of the appendix. By definition, the 
value of the maximum of /3,4[p, m] is the Legendre transformation of the functional logS[M] at point p, i.e. we 

have 



(iA[p] = sup PA[p, u] 



(3.7) 

= l3A[p,u], (3.8) 
where u is the solution of &A[p, u]/5u{x) = for a given p{x). If it exists, u is therefore the unique solution of 



p{x) 



5\og^[u] 



5u{x) 



(3.9) 



As usual, existence is much more difficult to prove than unicity. Clearly, if for some x e V, p{x) < Eq. (^^) will not 
have any solutions in U. It is easy to prove that, for the ideal gas {W^ = VA^) the condition p{x) > is sufficient 
for Eq. ( |3.9|) to have a solution. For non-trivial models, this condition is in general not sufficient since there exists 
an u pper bound on || p ||^i. However, the beautiful result that for H- stable systiams the set TZ of densities for which 
Eq. (3.9) have a solution is a convex set has been proved by Chayes and Chayes.Q 

Some important properties of the functional /3 ^[p] can easily be deduced from its mere definition. Firstly, note 
Young's inequality which follows readily from Eq. (BJ) 



(Vu e U) (Vp e TZ) \ogE[u] + l3A[p] > {p\u) 



(3.10) 



Secondly, f3A is a strictly convex functional of p. 
< A < 1. Since TZ is convex p\ = Api + (1 — 
corresponding to px- It follows from Eq. (pjh that 



To prove that, let pi,P2 7^ Pi two distinct densities of TZ and 
A)p2 is also a density of TZ. Let ma the solution of Eq. (|3.9|) 
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l3A[\pi + (1 - X)P2] = X(3A[ux,pi] + (1 - A)M[ma, P2] 
< \l3A[pi] + (1 - \)PA[p2] , 



(3.11) 



which proves the strict convexity. An alternative proof of the convexity of can be found in Ref.Ll Finally the 

functional derivative of PA[p\ with respect to p is easily obtained from Eq. ( |3.8| ). We have 

S/3A[p] SpA[p,u] 



6p{x) 5p{x) ' 

= w(^) + (P I 
= u{x) , 



Su ^ log S 
Sp{x) Su 



6p{x) 



where we have made use of Eq. (3.£) 



Note that it follows from Eqs. (2.6) and (3.8) that we can rewrite 



(3.12) 



(3.13) 



where /o(w) denotes the grand canonical equilibrium density in the phase space for the system in the presence of 
the external potential u.rjTChe expression (3.13) of /3^[p] coincides of course with that of the intrinsic Helmholtz 
(Kohn-Sham) free energyHiZl 
We would like now to comment briefly on Eqs. ( |3.7D , ( |3.8| ). They can be rewritten more explicitly as 

(Vu(f) ) f3A[p] > [ d'^x p{x)u{x) + I3n[u] = (3A[p,u] , (3.14) 



where the equality holds for the unique potential u = u solution of Eq. (3.9) . Eq. ( 3.14| ) takes thus the form of a 
variational principle which asserts that the generalized potential u{x) which gives rise to a given density p{x) is that 
which maximises the functional l3A[p,u]. If it exists it is unique. We have thus obtained a formal solution of the 
inverse problem to that usually considered in the DFT. To be effective this variational principle requires an exact or 
at least an approximate knowledge of the functioiafiL/3il[u]. Unfortunately, although many efforts have been devoted 
in the past to find good approximations of /J^fpJElLl for many fluids of interest, similar approximates for /3ri[u] are 
not available, at least to the author's knowledge. P-l 
The functional /3^[p] is the generating functional of the so-called direct correlation hmctions c''"-'(a;i, . . _a;„).U 
These functions play in the theory of liquids the role devoted to the vertex functions in statistical fleld theoryJlj They 
are defined as 



(xi, , 



1 •^ri. ) 



(5(") l3A[p\ 
5p{xi) ■ ■ ■ 5p{xn) 



(3.15) 



The direct correlation fuBe|tians c'"-* and the UrseU correlation functions ci""* are related through generalized 
Ornstein-Zernike relations.llirll3 In fact, for historical reasons, one rather defines the "true" direct correlation functions 



c(")(xi, . . . ,f„) by 



c(")(i?i, 



! •^n ) 



(fi 



(3.16) 



where the c|^'' are the c*^"^ functions of the ideal gas, so that the c'-"-' are the functional derivatives of minus the excess 
free energy. It is easily established that 



1 '^n) 



-log (A^pix)) 



- (n-2)! J](5'^(fi -f,) (n>2). 



(3.17) 



so that, for instance we have \xi, X2) — {xi, X2) + S {xi — X2)/ p{xi). Since the functional is convex, it 

(2) 

follows by application of the theorem II of the appendix that the quadratic form PAp [h] , where h in an arbitrary no 
zero function defined on V, is definite positive, which can be expressed as 



d''xid'^X2 



5'\xi 



X2 



p{xi 



^^Hxi,X2)]Hxi)h{x2)>0^ 



(3.18) 



Note that, in the above equation, a priori the inequality is non strict (see the appendix for more details on this point). 
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C. The DFT recovered 



As for functions of several variablesEB, the Legendre transform of a functional is involutive, i. e., in our case, 
the Legendre transform of /3^[p] coincides with logS[M]. Indeed, since is a strictly convex functional defined 

on 7?., one can define its Legendre transform /3A[it] for any u € U. One first defines the two fields functional 



(3.19) 



which is concave in p for a given u. The Legendre transform /3A[u] is then defined as the maximum of /3A[m, p] at u 
fixed, i.e. 



(3A[u] — sup f3A[u,p] , 
pen 



(3.20) 

= /3A[ii,p], (3.21) 
where p is solution of SA[u, p]/Sp(x) = for the given field u{x). If it exists, p is therefore the unique solution of 

S(3A[p] 



u{x) 



5p{x) 



(3.22) 



In order to prove that {3 A = log^ we insert the expression ( |3.8| ) of /9^[p] in the r.h.s of Eq. ( 3.21 ) which yields 

(3A[u] = {-p\{u - u)) + logS[u] , (3.23) 

where u is the unique solution of 

(5 logs 



p{x) 



5u{x) 



(3.24) 



We claim now that u = u and therefore that (3 A = log^. Indeed it follows from Eq. ( 3.12 ) that 

5pA[p\ 



u{x) 



5p{x) 



p=p 



and hence, by comparison with Eq. (3.22) u = u, Q ED. 



We would like to comment briefly on Eqs. (3.20), (3.21). They can be recast under the more familiar form 
(Vp(a;)) pn[u]<PA[p]^ [ d'^x p{x)u{x) = (3nv[u, p]{= -(3A[u, p]) , 



(3.25) 



where the equality holds for the unique profile p = p solution of Eq. ( 3.22 ). Eq. ( 3.25 ) is precisely the DFT 
variational principle which states that the density p(x) coixesponding to the imposed external potential u{x) is that 
which minimizes the functional (3ftv[p, u] for a given M(a;).El Q Moreover the solution is unique (if it exists). 



IV. SUMMARY AND ADDITIONAL COMMENTS 



In this paper we have studied some of the properties of a nonuniform classical fluid occupying a finite volume V 
and in the presence of an external potential (p{x). The thermodynamic and structural properties of the fluid are 
described by the logarithm of the grand partition function logS[u] which is the generating functional of the connected 
correlation fonctions g[^\ logS[M] is a convex functional of the generalized external potential u{x) = (3p — f3ip{x) 
and its Legendre transform /?^[p], which identifies with the Kohn-Sham free energy, is a convex functional of the 
density p{x). PA[p] is the generating functional of the direct correlation (or vertex) functions c'"^ The functions 
Gi"^ and c(") are related by generalized Ornstein-Zernike relations. The Legendre transform is involutive and thus 
the Legendre transform of /3^[p] coincides with logS[it]. As a consequence of the convexity of /9^[p] and logS[u] the 
two following quadratic forms are definite positive 

logSf[/i]>0, (4.1) 
PA^P[h]>0. (4.2) 
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Moreover, for iJ-stable systems, the set U of the generahzed external potentials u and the set TZ of the densities p are 
both convex sets. 

For arbirary potentials u d U and densities p £ TZ we have Young inequality 

^A[p]+\ogE[u]>{p\u). (4.3) 
/3„4[p] and logS[u] are both the unique solution of a variational principle : 

pA[p]^sup{{p\u)~logE[u]) (Vpe7^), (4.4) 

ueu 

\ogE[u] = sup {{p\u) - l3A[p]) (Vw e U) . (4.5) 

pen 



The latter variational principle (4^) is that of the DFT whereas the dual one (4^) is a new one. 

All the results which are summarized above are valid for any finite domain V, before the passage to the thermody- 
namic limit, and one can wonder whether they survive in the infinite volume limit. The notion of thermodynamic limit 
is not well defined for an inhomogeneous system and we concentrate henceforth on the case of homogeneous systems. 
The domain V that we consider now could be a cube of volume V — L"^ with periodic boundary conditions, or the 
surface of a hypersphere, etc. There is no applied external potential and u = Pp reduces to the chemical potential, so 
that our system is homogeneous even for a finite volume. The grand-canonical pressure pv is defined by the thermo- 
dynamic relation logS(/i) = f]pvV, and similarly we define ay(p) = A{p)/V. Of course all the results of the paper 
apply to this pecular case and one concludes that, for a finite volume V, the pressure pvifJ-) at given temperature and 
volume is a strictly convex function of the chemical potential and that the specific grand-canonical free energy ay (p) 
is a strictly convex function of the density. This is true even in the event of a liquid-vapor transition for instance. 
This behavior of av{p) contrasts with that of the canonical free energy a™ "(p) which is a bimodal function of p in 
the two phases region, due to surface effects. But, what happens in the thermodynamic limit? We first note that the 
strict inequalities which enter the definition of strict convexity for both pv{l^) and ay(p) can become equalities in 
the infinite volume limit. Therefore, in the TL, we can only claim that, at a given temperature, Poc{p) is a convex 
function of p and that aoo{p) is a convex function of the density. Moreover, as it is well known, the continuity of 
some derivatives can also be lost in the TL. This is the case for Poo{p) which exhibits a "cusp" for some p{T) at a 
temperature T below the critical temperature Tc- One can however still define a Legendre transform in case of such 
discontinuities. Let us denote by Pg and pi the two distinct values of dpoo/dp at the cusp, then Ooq^plwill be a linear 
function of p in the interval [pg, pi]. Cusps and intervalls correspond in a Legendre transformation.ErOJ To summarize, 
in the coexistence region, the strict convexity of av{p) is lost in the TL ; by contrast, poo{p) remains strictly convex 
in the TL even in the event of a phase transition, but its derivatives have discontinuity points. 
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APPENDIX A: CONVEX FUNCTIONALS 

Functionals are a simple generalization of functions of several variables; they can be-. defined as "functions of 
functions" . In this appendix we give formal proofs of the theorems on convex functionalstj which are needed in the 
text, without aiming at any mathematical rigor. The ppapfs given below are a mere transcription of the material of 
any standard textbook on functions of several variables£3 

We consider a convex set X of real functions defined on a domain V C M'' with a norm || Lp \\ {tp £ X). We define 
a functional -F as a mathematical object which associates to each element ip a real number F[lp\. F will be said 
differentiable at point 1^9 of X if it exits a linear functional F^jP such that 

Fy + M = py] + [M + e[M II h II , (Ai) 

lim e[M=0, (A2) 



where we have noted 



SF 

F«[M= / rf'^^T^M^), (A3) 



V 



6lp{x) 
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where 5F/5(p{x) is the functional derivative of F with respect to the field (p{x). One deduces from Eq. ( Al) a useful 
expression for the derivative of F in the direction 6ip : 

Fi^^ m ^ lim^ ^[^ + AM-^M ^ ^) ^ 

A — '0 A 

It is easy to prove by applying RoUe's theorem to the real function g{t) = F{lp + tSip) that aTaylor Mac-Laurin 
expansion also holds for functionals : 

F[^ + S^]^F[ip]+Fl%^M (0<A<1). (A5) 
If F is twice diffcrentiable it generalizes as 

F[^ + S^] = F[^]+ [S^] + If^^I.s^ M (0 < a < 1) , (A6) 

where 

It will be necessary to distinguish carefully between and strictly convex functionals. F[ip\ will be said convex if 
(V(^i, ip2 G X) and for all < A < 1 one has F[Xipi + (1 — X)'P2] < ^F[fi] + (1 ~ ^)F['P2], whereas it will be said 
strictly convex if (V</5i, <p2 ^ (fi <^ X) and for aU < A < 1 one has F[Xipi + (1 - A)iy92] < XF[ipi] + (1 — X)F[ip2]- F[ip] 
will be said concave (resp. strictly concave) if —F[ip] is convex (resp. strictly convex). Henceforth we shall assume 
moreover that all the considered functionals are at least twice differentiable. Let us establish our first theorem which 
states that a convex functional lies above its tangent plane, i.e. more precisely 
Theorem I : 

• F 

^ (V^i, </^2 G X) F[ip2] > F[ipi] + F^\\ip2 - ipi] . 

• F strictly convex <^=> (V(y52 =/= fi E X) F[ip2] > F[ipi] + F^]^[ip2 - tpi] . 

Suppose F convex and let < A < 1 and ((pi, (p2 7^ ^1 ^ X). The convexity of F implies that 



F[ipi + \{ip2 - (pi)] - F[ip 



- < F[^2] - F[^i] 



A 



Taking the limit A ^ and making use of Eq. (A4) yields F[(/32] > -^'[(pi] + F!^-^ [ip2 — ifii]- If F is strictly convex the 
above proof does not imply the strict inequality, more care is needed. Let < lu < 1. Since 

(1 - X)(pi + X(f2 = (Pl + — i0{ip2 - (Pl)) , 

for < A < w we have 0<((jJ — A)/ijj<1 and therefore 

F[lpi + A((^2 - ^1)] < - — -F[ipi] + —F[ipi + Lo{ip2 - ipi)] ■ 

If F is strictly convex one thus have for < A < 

F[lpi ^ X{ip2 - ipi)] - F[ipi] F[lpi + uj{if2 - Vi)] - F[lpi] 

< < F[ip2\ - F[ipi\ . 

A U! 

Taking the limit A ^ and making use of Eq. ( A4) yields the strict inequality : F[ip2] > F[ipi] + F^j' [ip2 — tpi], since, 
by hypothesis, a; < 1. 

Reciprocally, we assume that 

(V./71, (^2 G X) F[ip2] > F[^i] + F«[(p2 - ^1] ■ 
Let ipi and ip2 two distinct points of X and < A < 1, we have 
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F^i] > F[X^, + (1 - A)^2] + i^S+(i-A)^. [(1 - - V2)] , (A8) 

F[^2] > F[A^i + (1 - A)^2] + [A(^2 " ^i)] . (A9) 

It is now sufficient to add the two above inequalities, multiplied respectively by A and by 1 — A to obtain 

XF[^i] + (1 - A)F[(^2] > F[\^i + (1 - A)(p2] , 

which establishes the convexity of F or its strict convexity if the inequalities are strict. 
Let us prove now the Theorem II : 

(2) 

• F convex 4=^ (V(y9i, if2 e X) F^^'[(p2 -'fi]>0, 

(2) 

• (V^?!, </32 ^ G X) Fi^^'[(p2 — fi] > =4> F strictly convex. 

Let us first assume that F is convex and let (pi an arbitrary function of X. We define G[lp] — F[ip] — F!^Wip — ipi]. It 
follows from the convexity of F and from theorem I that {Vip G X) G[(p] — G[(pi] > 0. We perform now a 2nd order 
Taylor Mac-Laurin expansion of G about (pi (cf. Eq. (A6) ) which yields 

G[p^+tp] G[p] = I II ^ > , 

lime(t) = 0. (AlO) 

Taking the limit t yields F^f [ip] > 0. 

Reciprocally, let ipi and ip2 {^2 7^ Pi) two arbitrary functions of X. It must exist some p g]</3i,</32[ such that (cf. 
Eq. (A6)) 

F[p2\ = F[p^] + [p2 - Pi] + [p2 - Pi] ■ (All) 



We define the (strictly) positive number by (^2 ^ = pi'^ ^ ^i) which ahows us to rewrite (All) as 

2 

F[p2] - F[pi\ - F^^^[p2 - Pi] - y4'^[^ - ^i] • 

Assuming the positiveness of the quadratic form of the r.h.s. yields, by application of theorem I, to the convexity of 
F and to its strict convexity in the case of strict positiveness. 
The last result we need is Theorem III : 

• If _F is convex and has a local minimum aX pi ^ X , then the minimum is global. 

• If is strictly convex then F has at most one minimum and this is a strict minimum. 
Let ^p — pi ~\- 5p an arbitrary function in X and < A < 1. Since F is convex 

F[pi + A^^] < (1 - X)F[p>i] + \F[p] ■ 

hence F[ipi + \5p] — F[pi] < X{F[p] — F[pi]). Assuming the existence of a local minimum of F at ipi implies that 
for some < Aq < 1 we have FYpi + X^Sp] — F[pi] > from which it follows that F[p] — F[pi] > 0. Therefore the 
minimum is global. If F is strictly convex the same demonstration leads to the inequalities 

< F[pi + XoSp] - F[pi] < Ao {F[p] - F[pi]) , 

which prove that the minimum is strict and therefore unique. 
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